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FOREWORD 


Thi«  report  presents  work  which  was  performed  under  the  Joint  Army  Navy 
Aircraft  Instrumentation  Research  (JANAIR)  Program,  a  research  and  exploratory 
development  program  directed  by  the  United  States  Navy,  Office  of  Naval 
Research.  Special  guidance  is  provided  to  the  program  for  the  Army  Elec¬ 
tronics  Command,  the  Naval  Air  Systems  Command,  and  the  Office  of  Naval 
Research  through  an  organization  known  as  the  JANAIR  Working  Group.  The 
Working  Group  is  currently  composed  of  representatives  from  the  following 
Offices: 


U.  S.  Navy,  Office  of  Naval  Research 
Aeronautics,  Code  46l,  Washington,  D.  C. 

-  Aircraft  Instrumentation  and  Control  Program  Area 

U.  S.  Navy,  Naval  Air  Systems  Command 
Washington,  D.  C. 

-  Avionics  Division;  Navigation  Instrumentation  and 
Display  Branch  (NAVAIR  5337) 

-  Crew  Systems  Division;  Cockpit/Cabin  Requirements 
and  Standards  Branch  (NAVAIR  5313) 

U.  S.  Army,  Army  Electronics  Command 
Avionics  Laboratory,  Fort  Monmouth,  New  Jersey 

-  Instrumentation  Technical  Area  (AMSEL-VL-l) 


The  Joint  Army  Navy  Aircraft  Instrumentation  Research  Program  objective  is: 
To  conduct  applied  research  using  analytical  and  experimental  investigations 
for  identifying,  defining  and  validating  advanced  concepts  which  may  be 
applied  to  future,  improved  Naval  and  Army  aircraft  instrumentation  systems. 
This  includes  sensing  elements,  data  processors,  displays,  controls  and 
man/machine  interfaces  for  fixed  and  rotary  wing  aircraft  for  all  flight 
regimes . 
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NOTICE 


Change  of  Address 

Organizations  receiving  JANAIR  Reports  on  the  initial  distribution  list 
should  confirm  correct  address.  This  list  is  located  at  the  end  of  the 
report  just  prior  to  the  DDC  Form  1473.  Any  change  in  address  or 
distribution  list  should  be  conveyed  to  the  Office  of  Naval  Research, 
Code  46l,  Washington,  D.  C.  20360,  Attn:  JANAIR  Chairman. 

Disposition 

When  this  report  is  no  longer  needed,  it  may  be  transmitted  to  other 
organizations.  Do  not  return  it  to  the  originator  or  the  monitoring 
office . 

Disclaimer 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  Defense  or  Military  Department  position  unless  so 
designated  by  other  official  documents. 
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ABSTRACT 


A  mathematical  model  that  positions  and  moves  a  variable  sised  23  pin 
joint  articulated  stick-man  in  a  crewstation  environment  is  presented. 
The  model  simulates  the  motion  of  pilots  in  a  given  cockpit  configura¬ 
tion  after  testing  the  gross  reach  capability  required  by  a  task.  It 
utilizes  a  non-linear  optimization  technique  to  position  and  orient 
the  joints,  analyzes  the  viewing  capability  after  the  operation  and 
detects  body  intersections  with  the  seatback  during  the  task. 
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1.0 


INTRODUCTION 


The  Cockpit  Geometry  Evaluation  project  has  as  its  goal  the 
development  of  a  digital  computer  program  which  would  evaluate 
the  geometry  of  proposed  crewstation  designs.  The  man-model 
will  be  used  to  evaluate  the  geometric  configuration  of  the 
crewstation  in  terms  of  the  crew's  physical  characteristics 
and  capabilities. 

PHASE  REQUIREMENTS 

There  will  be  six  phases  to  the  development.  Each,  of  one  year's 
duration,  will  result  in  a  successively  greater  detailed  descrip¬ 
tion  of  anthropometric  characteristics  and  capabilities  of  the 
pilot  (Figure  l) . 

In  Phase  I,  a  "baseline"  man-model  will  be  developed,  using  a 
skeletal  frame  consisting  of  stick  hr  line  links  intersecting 
at  pin" -joints  with  specified  kinds  of  link  movement 
originating  from  the  joints. 

To  simplify  the  problem  for  Phase  I,  attention  is  restricted  to 
an  aircraft  pilot's  station.  For  the  present  development,  the 
copilot  will  be  immobile  and  of  fixed  dimension.  Furthermore, 
both  pilot  and  copilot  will  be  operating  from  a  seated  position. 

In  Phase  II,  skin  volumes  of  rectilinear  or  circular  section 
will  be  added  to  the  above  "pin-joint  stick-man." 

Phase  III  will  provide  some  force  and  control  grasping  capability 
to  the  "pilot"  and  incorporate  the  "preferred"  positions  for 
control  grasping  and  force  application. 

Phase  IV  will  allow  for  finger  dexterity  and  excursion 
capability  to  specified  joints  of  the  pilot. 
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Figure  1.  6-Year  Plan 
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In  Ihase  V,  additional  force  vector  variables  and  initial 
energy  expenditure  variables  will  be  incorporated  into  the 
"pilot's"  range  of  capability. 

Finally,  Phase  VI  will  provide  a  "skin"  deformation  capability 
and  the  final  energy  expenditure  variables  will  be  added  to 
the  model. 

During  each  of  the  Phases,  refinements  in  definitions  of  cockpit 
dimensions,  display-control  locations  and  control  shapes  will 
be  made.  In  addition,  validations  will  be  performed  to  test  the 
hypothesis  that  the  pilot  model's  positions  and  movements  do 
not  depart  "significantly"  from  those  of  human  pilots. 

EVALUATION  DESCRIPTION 

In  order  to  make  the  evaluation  more  meaningful,  the  computerized 
man-model  (when  all  project  phases  are  complete)  should  be  able 
to  evaluate  a  cockpit  design  in  terms  of: 

(1)  Any  human  anthropometric  combination  of  link  sizes  between 
the  first  and  ninety-ninth  percentile 

(2)  Any  set  of  pilot  task  sequences 

The  computerized  man-model  will  accommodate  any  set  of  cockpit 
dimensions  and  control  shapes  as  well  as  any  set  of  cockpit- 
induced  new  movement  restrictions  due  to  seats,  harnesses, 
panels  and  clothing. 

In  addition,  the  computerized  man-model  will  supply  constructive, 
consistent,  and  reproducible  predictions  concerning: 

(1)  Display  and  Control  Accessibility  -  Is  it  possible  for  the 
"pilot"  to  view  a  cockpit  display  (such  as  an  altimeter  on 
the  control  panel  or  a  specific  area  in  the  windscreen) 
and  can  he  also  "see"  what  he  is  doing  while  reaching  for 
a  control? 


D162-10128-1 

3 


(2)  Control  Actuation  -  Is  it  possible  for  the  "pilot"  to 
operate  or  can  his  appendages  be  placed  in  the  proper 
orientation  for  "handling"  the  control,  given  the 
accessibility  of  the  control? 

(3)  Pilot  Workload  -  Over  any  given  set  of  tasks,  how  much 
"energy"  is  expended  by  the  pilot;  hov  much  "work"  is  done 
by  the  pilot  in  terms  of  the  body  mass  displacements  and 
distances  traveled  by  each  link;  and  what  is  the  amount 

of  visual  "activity"  as  measured  by  head  and  eye  deflections? 

(4)  Interference  -  Is  there  physical  or  visual  interference 
occurring  between  the  pilot  and  the  cockpit  structure  or 
equipment  during  a  task  to  be  performed?  The  pilot's 
body  and  personal  equipment  will  also  be  considered  as 
possible  obstructions. 

(5)  Suggested  Areas  for  Design  Change  -  In  order  to  make  the 
cockpit  physically  compatible  with  pilot  capabilities,  the 
computerized  man-model  will  specify  why  a  task  is  infeasible. 
This  will  suggest  that  a  change  be  made  in  the  cockpit 
design  or  procedures. 
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2.0 


SUMMARY 


2.1  PROJECT  DESCRIPTION 

The  problem  is  to  develop,  over  a  period  of  six  years,  a 
computerized  man-model  as  a  means  of  evaluating  competing  designs 
of  aircraft  crewstations .  Specifically,  this  is  to  be  accomplished 
by  determining  the  ability  of  any  size  of  crewman  to  reach  or 
view  the  controls  and  displays  of  his  crewstatlon.  Interference 
with  projecting  parts  of  the  cockpit  during  the  motion  is  also 
to  be  determined.  These  determinations  are  to  be  made  by  simulating 
the  movements  of  a  human  pilot  during  the  execution  of  various 
tasks  in  the  crewstatlon  environment  by  means  of  a  mathematical 
model  programmed  for  the  computer.  When  given  the  trajectories 
and/or  orientations  of  the  pilot's  extremities  (e.g.,  the  hands), 
this  model  is  to  predict  the  trajectories  and  orientations  of 
the  remaining  body  segments. 

The  project  objectives  for  Phase  I,  therefore,  are  the  following: 

a)  Develop  a  math-model  and  an  associated  computer  program  to 
position  and  move  a  23-pin- joint  stick-man  in  three  dimensions. 

b)  Develop  validation  procedures  and  criteria  to  assess  the 
’reality"  of  model  positions  and  movements  compared  to 
actual  pilot  movements. 

c)  Satisfy  the  requirements  of  Phase  I  of  the  project  by 
supplying  useful  predictions  concerning: 

o  Display-control  accessibility 
o  Control  actuation 

o  Workload  in  terms  of  energy  expended,  body  mass 
displacements,  and  distances  traveled  by  the 
flight  crew's  eyes  and  each  flight  crew  appendage 
o  Visual  interferences  between  the  flight  crew  and 
the  cockpit  structure /equipment  and  physical 
interferences  between  flight  crew  appendages  and 
the  seat back. 
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o  Design  changes  required  to  make  the  cockpit  subsystem 
compatible  with  flight-crew  physical  capabilities. 

RESULTS 

The  present  model  relies  on  optimization  theory  to  find  the  body 
segment  trajectories  and  orientations  by  minimizing  a  cost 
(objective)  function.  This  cost  function  is  assumed  to  express 
the  amount  of  effort  required  by  the  crew-man  to  achieve  a  given 
body  configuration.  Total  mass  displacement  and  visual  inter¬ 
ference  that  occur  during  the  task  motion  are  analyzed  after  the 
simulation  of  a  task. 

A  version  of  the  man  motion-model  that  simulates  a  human  torso 
and  right  arm  system  has  been  completed.  An  experimental  motion 
model  for  the  upper  body  (torso,  arms,  and  head)  has  also  been 
developed.  The  current  upper  body  model  provides  for  all  body 
structural  constraints  except  the  joint  angular  limits. 

CONCLUSIONS 

The  end  result  of  the  Phase  I  Cockpit  Geometry  effort  is  a  work¬ 
ing  experimental  model  of  human  motion.  It  appears  suitable 
for  use  as  a  comparative  tool  in  limited  cockpit  geometry 
evaluation  studies.  It  also  appears  suitable  as  a  baseline  for 
further  research  and  development  aimed  at  obtaining  the  fully 
applicable  man-model  desired. 
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2.4 


RECOMMENDATIONS 


(1)  The  design  of  the  present  model  reflects  what  may  be 
a  highly  simplified  view  of  the  way  human  pilots  move 
during  execution  of  a  task.  Therefore  additional  data  on 
human  motion  should  be  gathered  to  aid  in  the  formulation 
of  better  mathematical  and  physical  hypotheses  concerning 
human  motion. 

(2)  More  research  on  the  modeling  problem  is  needed  to  Insure 
accuracy  and  stability  in  the  model.  The  baseline  man- 
model  is  flexible  enough  to  incorporate  changes  in 
mathematical  or  physical  hypotheses  that  might  result  from 
this  research.  Pull  advantage  of  this  flexibility  should 
be  taken  by  devoting  much  of  the  Phase  II  effort  toward 
major  improvements  in  motion-model  formulation. 
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3.0 


DISCUSSION 


3.1  ASSUMPTIONS  AND  NOTATION 

The  labelling  of  points  on  the  stick-man  is  presented  in  Figure 
2  for  the  purposes  of  this  document. 

Task  motions  are  performed  in  an  incremental  fashion,  a  task 
being  defined  by  moving  terminal  body  segments  (e.g.,  the 
hands)  a  step  at  a  time  toward  the  control  points  for  these 
segments.  Solutions  are  obtained  for  the  Euler  angles  that 
define  the  stick-man's  configuration  at  each  step  by  using 
optimization  with  non-linear  constraints.  Trajectories  for 
the  intermediate  body  segments  at  each  step  are  calculated  from 
the  Euler  angles.  Here,  an  underlying  assumption  is  that  human 
motion  can  be  simulated  by  obtaining  solutions  in  a  step-wise 
manner  and  using  interpolation  to  obtain  body  configurations 
intermediate  to  the  calculated  ones.  Also,  since  time  is  not 
actually  a  parameter  in  the  present  model,  it  is  assumed  that  the 
motion  is  performed  with  constant  velocity  so  that  discrete 
time  values  may  be  assigned  to  positions  along  the  motion 
trajectory  of  the  body. 

The  main  assumptions  of  the  overall  model  and  the  notation  that 
applies  to  Figure  2  are  as  follows: 

ASSUMPTIONS 

(l)  The  stick-man  (Figure  2)  is  initially  placed  in  the 

cockpit  so  that  the  eye  reference  point  coincides  with 
location  P^.  Since  initial  link  lengths  and  orientations 
are  given  (with  respect  to  the  preceding  Joint  or  point), 

the  other  P.  are  calculable  with  respect  to  P  . 

J  7 
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Figure  2.  23  Joint  Man-Model  Link  Nomenclature 


(2)  In  particular,  the  location  P  is  found  and  P  becomes 

o  o 

the  fixed  origin  for  all  movements  of  the  stick-man 
thereafter.  (Thus  the  locations  Pj  are  again  determined 
with  respect  to  the  new  origin.  All  cockpit  locations 
are  calculated  with  respect  to  P^.) 

(3)  A  task  is  given  (i.e.,  a  specific  total  body  movement 
from  an  initial  terminal  joint  position  to  a  final 
terminal  joint  position)  specifying: 

a)  Final  terminal  joint  locations  (each  P^*). 

b)  New  orientations  for  each  terminal  link  of  length  L^*. 

c)  A  time  value  to  achieve  the  final  position. 

d)  A  time  value  to  remain  in  the  final  position. 

e)  A  cockpit  location  or  plane  area  at  which  the  eyes 
are  to  focus. 

(4)  The  legs  can  be  ignored  for  the  Phase  I  model,  leaving 
only  the  upper  torso,  head,  and  arms  to  consider. 

NOTATION 

(1)  Let  P  denote  the  location  relative  to  P  of  point  j,  where 

v  O 

j  ranges  from  j  =  0  to  j  =  34  (e.g.,  Pq  =  (0,0,0)). 

(2)  Let  L.^  denote  the  length  of  Link  i.  Link  i  is  the  line  or 
polygonal  line  connector  between  two  adjacent  junctions 
(points  and/or  joints). 

(3)  For  terminal  points,  j  is  restricted  to  the  values 

j  -  6,  7,  8,  21,  22,  31,  33,  32,  34  only  (for  j  -  6, 

7,  8,  21,  31,  33,  32,  34,  j  is  an  interior  point). 

(4)  Let  Pg  denote  the  location  relative  to  P^  at  the  top  of 
spine  point  defined  as  the  midpoint  of  the  line  (P^,P10) 
intersecting  the  polygonal  line  (P^P^). 

(5)  Let  P^  denote  the  location  of  the  junction  occurring 
between  joints  Pg  and  P^. 


D162-10128-1 

13 


(6)  Let  Pj.  denote  the  location  of  the  junction  occurring  at 
the  intersection  of  the  polygonal  lines  (P^,Pg)  and 

(7)  Let  P^  denote  the  midpoint  of  the  line  (Pg,Pg)  between 
the  eye  joints. 

(8)  Let  (Pg2,P^)  be  the  line  with  length  L^  +  L^  and 
endpoints  P^  (toe) ,  P^2  (heel) ,  and  with  midpoint  P^ 
(right  foot); 

Let  (*3! >*33)  ^  ^he  line  with  length  L^  +  L^  and 
endpoints  P^^  (toe) ,  P^  (heel) ,  and  with  midpoint 
(left  foot). 

REACH  ANALYSIS 

The  Reach  Analysis  is  concerned  with  determining,  in  a  gross 
manner,  if  a  prescribed  task  is  with  BOEMAN-I's  (Figure  3) 
reach.  That  is,  are  BOEMAN-I's  link  dimensions  and  joint 
angular  limits  such  that  he  can  simultaneously  position  the 
palm  joints  of  each  hand  at  the  task  prescribed  locations 
without  rising  from  the  seated  position,  (it  is  assumed  that 
the  lumbar  joint  (JD)  is  fixed  and  immobile.) 

To  this  end,  two  problems  are  formulated: 

(a)  Is  the  task  feasible? 

(b)  If  not,  redefine  the  task  control  locations 
for  the  hands  to  insure  task  feasibility. 

The  solution  to  the  latter  problem  allows  the  motion  model  to 
proceed  and  identify  any  additional  difficulties  that  could 
occur  such  as  in  body  orientation  or  viewing  ability. 
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Figure  3.  23-Joint  Man-Mode/  Link  Nomenclature 
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3.2.1  Task  Feasibility 

It  is  natural  to  consider  the  top  of  the  spine  (J^)  as  the 
key  to  task  feasibility,  for  it  is  a  juncture  of  three  body 
systems:  the  spine,  the  right  arm,  and  the  left  arm.  By 
fixing  the  palm  joints  at  the  task  locations,  with  the  lumbar 
joint  already  fixed  relative  to  the  seat  reference  point,  the 
task  is  deemed  feasible  if  these  body  systems  yield  at  least 
one  common  top  of  spine  location. 

Consequently,  three  regions  will  be  defined  with  respect  to 
the  fixed  points  of  each  body  system,  each  describing  the 
locus  of  points  that  the  top  of  the  spine  can  assume.  Then 
task  infeasibility  will  be  assumed  if  it  is  not  possible  to 
exhibit  a  point  in  the  mutual  intersection  of  the  three  regions. 

HAND  REGIONS 

The  Right  Hand  Region,  denoted  Rl,  is  defined  as  the  set  of 
all  points  that  the  top  of  the  spine  can  assume,  given  that 
the  right  palm  joint  is  at  the  corresponding  task  control 
point.  It  is  defined  by  the  two  extreme  arm  positions,  encom¬ 
passes  virtually  all  possible  arm  positions  (see  Figure  5). 

The  left  hand  region,  denoted  R2,  is  defined  in  a  similar  manner 
(see  Figure  6). 

For  BOEMAN-I,  given  a  fixed  palm  location,  the  effective 
distance  between  the  palm  and  any  top  of  spine  location  ranges 
from  zero  to  the  sum  of  the  link  lengths  for  the  corresponding 
arm  system  (including  the  shoulder) .  It  is  necessary  to 
determine  bounds  on  this  distance  that  closely  approximate 
those  encountered  by  human  beings. 
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For  a  lower  bound,  it  is  assumed  that  some  positive  radius, 
r1,  equal  to  one  half  the  length  of  the  interclaricular  link 
(about  an  inch)  be  used,  where  i  -  1,  2  denote  the  left  and 
right  am  system,  respectively. 

(1)  r±  «  ,  i  =  1,2 

A  small  positive  distance  is  assumed  because  it  is  not  possible 
to  touch  the  top  of  the  spine  with  either  palm  due  to  the 
intervening  flesh. 

For  an  upper  bound,  it  is  necessary  to  choose  a  radius  that 
represents  a  compromise  between  forward  am  reach  lateral 
am  reach.  In  the  former,  the  links  between  the  top  of  the 
spine  and  the  shoulder  do  not  contribute  to  the  for¬ 

ward  reach;  however,  they  contribute  fully  in  the  latter. 

This  holds  of  course,  when  one  aeglects  the  rotational  capability 
of  the  spine.  However,  a  twisting  spine  (torso)  tends  to  give 
added  reach  to  only  one  am  system  and  compensates  by  decreas¬ 
ing  the  reach  of  the  other  am  system.  Thus  to  approximate 
maximum  reach  of  human  beings,  and  effect  a  "balanced"  reach 

that  the  maximum  distance  is  given 

+  *13  +  *15  +  *17  +  *19)  ) 

*  *14  *  *16  +  ‘18  +  *20 )  / 

(See  Figure  4  ). 


capability,  it  is  assumed 
by  radius  R. , 
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Figure  4.  Maximum  Reach  for  Each  Arm  System  (Top  View) 
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Thus ,  R^  corresponds  to  an  extended  arm  position  and  r^ 
corresponds  to  a  contracted  arm  position.  This  will  include 
virtually  all  other  arm  positions  since  they  must  have  a 
reach  distance  (from  the  palm  to  the  top  of  the  spine) 
between  and  r^ . 

The  right  hand  region  (Figure  5)  then,  can  be  described  as  the 
volume  between  the  two  concentric  spheres ,  both  centered  at 
the  right  palm  control  point,  (x?,y2,z2)  and  with  radii 
r? ,R2 .  Then 

(3)  R1  =  {(x,y,z)  |  r2  s  (x-x?)2  +  (y-y2)2  ♦  (i-ZjfjR^} 

Similarly,  the  left  hand  region  (Figure  6)  can  be  described  as 
the  volume  between  the  two  concentric  spheres,  both  centered 
at  the  left  palm  control  point  (x^y^z^  and  with  radii  r^,R^ 
Then 

(M  R2  =  {(x,y,z)  |  r2  s  (x-x^2  +  (y-yx)2  +  (z-z^R2} 

SPINE  REGION 

The  spine  region,  denoted  R3,  is  defined  as  the  set  of  all 
points  that  the  top  of  the  spine  can  assume  given  that  the 
lumbar  joint  (Jq)  is  fixed  and  that  the  links  in  the  spine 
system  may  not  violate  any  joint  angular  limits. 

R3  differs  from  R1  and  R2  in  that  it  is  a  more  accurate  (and 
more  complex)  region.  This  is  because  it  is  necessary  for 
this  region  to  contain  no  points  that  cannot  be  attained  by 
appropriate  rotation  angle  values  on  the  spine  joints,  with 
specified  link  lengths.  It  is  generated  la  Wm 
distinct  bounding  surfaces  C^,C2,C^. 
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Figure  5. 


Right-Hand  Region ,  R  1  —Cross  Section  of  a  Portion  of  Two  Concentric  Spheres 
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MAXIMUM  REACH 


LEFT-HAND 
CONTROL  POINT 

(xl'  21> 


TOP  OF 
SPINE 


Figure  6.  Left-Hand  Region,  R 2 -Cross  Section  of  a  Portion  of  Two  Concentric  Spheres 
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The  first  bounding  surface  (c^)  is  developed  by  considering 
the  effect  of  angular  deviation  in  the  lumber  Joint  Jq 
(Figure  7  )  with  zero  deviation  in  the  thoracic  joint  (j^) 


Let  R  =  A+B, 

the  sum  of  link  lengths  in  the  spine 

a  = 

maximum  angular  deviation  from  the  standard 
(upright)  position  (left  and  right) 

J2= 

(x,  y,  z) 

V 

(0,  0,0) 

Then  the  top  of  the  spine  point  is  on  the  surface  if  and 
only  if 

(5)  x2  +  y2  +  z2  =  R2 


Subject  to: 

-R  sin  a  &  x  §  R  sin  a 

Osys  R  sin  a 

R  cos  a  %  z  s  R 

Next,  consider  the  fixed  extreme  limits  that  the  middle 
Joint  ( J^)  can  deviate  and  the  angular  deviation  of  the 
bottom  joint  ( Jq) .  The  result  is  the  bounding  surface 
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Figure  7  Joint  Angular  Motion  at  Bottom  Joint— Upper  Link  at  Zero  Deviation 
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Figure  8.  Joint  Angular  Motion  at  Bottom  Joint  —  Upper  Link  in  Extreme  Deviation 
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Let  0  « 


f  = 


g  = 


r 


2 

r  = 


2 

r 


maximum  angular  deviation  from  the  standard 
(upright)  position  towards  x-axis 

B  sin  0 

B  cos  0 

radius  of  inner  curve 
f2  +  (g+A)2 

B2  sin2  0  +  (A  +  B  cos  0  )2 

sin 20  +  A2  +  2AB  cos  0  +  B2  cos 20 

2  2 

A  +  B  +  2AB  cos  0 


Then  the  top  of  the  spine  is  on  surface  if  only  if 


(6) 

subject  to 


2  2  2  2 
x  +  y  +  z  *  r 


-  r  cos  6sx?r  cos  6 
0  <y  sr  cos  6 
h  %  z  s  r 


where  h  =  r  sin  5  ,  6  =  90  - « -  V  , 

V  =  cos'1  ^  *-+  B  cosjj^ 

(using  the  law  of  cosines,  B2  «  A2  +  r2  -2Ar  cosy 

=  A2  +  (A2  +  2AB  cos  0  +  B2)  -2Ar  cosV 
r  cosy  =  A  +  B  cos  0  ) 


Finally,  consider  the  extreme  limits  that  the  bottom  joint  (j  ) 
can  deviate  and  the  angular  deviation  of  the  middle  Joint 
(j^).  The  result  is  boundary  surface  . 
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CROSS  SECTION  OF 
SURFACEl  C3 
(LEFT  SIDE) 


Figure  9.  Joint  Angular  Motion  at  Middle  Joint-Lower  Link  in  Extreme  Deviation 
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On  the  positive  x-axis,  the  top  of  the  spine  point  is  on 
surface  Cg  if 

(7)  (x-3%)2  +  (y-y^)2  +  (z-z^)2  *  B2 

where 

0  «  £  «  7T 

0  §  os  jt/2 

and  subject  to: 

R  sin  o  %  x  s  cos  <5 
R  sin  o;  «  ys  cos  6 
h  <  z  $  R  cos  a 


x4 

y4 

zk 


A  sin  a  cos  e 
A  sin  a  sine 
A  cos  a 


On  the  negative  x-axis,  the  top  of  the  spine  point  is  on 
surface  if 

(8)  (x-x^)2  +  (y-y3)2  +  (z-z3)2  -  B2 

where 

*  -A  sin  a  cos  e  |  0  %  e  %  tt 

y^  =  +A  sin  a  sin  e  j  0?as  x/2 

=  A  cos  a  J 

and  subject  to: 

-r  cos  6<xs  -R  sin  a 
R  sin  a  s  y  5  r  cos  6 

h  %  z  %  R  cos  01 

where  h  =  r  sin  6  (see  (  6  )) 

Then  Cg  and  together,  describe  the  spine  region.  It  has 
been  assumed  that  the  seated  pilot  is  constrained  by  a  seat  back. 
Thus  C1  are  defined  for  y?0  only.  Region  R^  is  depicted  in 
Figure  10. 
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Figure  10.  Region  3-Feasible  Region  for  Top-of-Spine  Locations 
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SPINE  REGION  PARTITION 

In  order  to  exhibit  a  point  in  the  Mutual  intersection,  it  is 
convenient  to  partition  R3  as  in  Figure  n  ,  into  two  Mutually 
exclusive  portions  P1,B2  such  that: 


Portion  PI  is  constrained  by 


1  r  <  x  -ty  +z  <  R  ,  , 

r  <  x  +y  +z  <  R 

1 

-R  siho’SXSR  sia al  \ 

j 

j(*>y>z) 

*  J(x,y*z) 

( 

| 

0  s  y  S  R  sin  Of  1  J 

0  s  y  ( 

R  cos  a  %  z  §  R  ' 

R  cosa  <  z  ' 

Portion  P2  is  constrained  by 


P2  =  (x,y,z) 


2  _  2  2 .  2 
v2 


r  <  x  +y“+z~  (x^  =  A  sin  a  cos? 


(x-x6)^  +  (y-y6)2  +  (z-z6)2<B2;  -  A  sinar  sinY 

-R  cos  5  R  cos  6  z^  ■  A  cos  a 

0  ?  ys  R  cos  <5  0  *ot<ir/2 

r  sin  6  «  zs  R  cos  a  0  s  y  I  ir  ) 


These  portions  allow  for  a  method  of  solution  using  an 
optimization  technique  (MINUM)  for  which  it  is  necessary  that 
upper  and  lower  bounds  be  available  for  each  of  the  unknown 
quantities  (x,y,z).  Hence,  the  overdetermined  constraint 
set  is  used  for  PI. 
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Figure  11.  Region  3  Partition 
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METHOD  OP  SOLUTION 


The  regions,  PI,  R2,  P3,  are  tested  for  a  mutual  intersection 
under  the  assumption  that  if  such  an  intersection  exists  a 
feasible  top  of  spine  location  is  calculable,  and  thus  the  task 
is  feasible. 

It  is  sufficient  to  show  that  the  mutual  intersection  region 
exists  by  calculating  a  point  belonging  to  all  three  regions 
simultaneously. 

Finding  a  feasible  top  of  spine  location  requires  that  there 

exists  a  point  (x,y,z)  in  Regions  1  and  2  and  in  at  least  one 

of  the  portions  PI,  P2  of  Region  3.  If  the  mutually  intersecting 

regions  have  a  finite  volume,  there  are  an  infinite  number  of 

feasible  candidates .  Thus,  it  is  necessary  to  select  one 

point  from  the  infinite  set.  To  do  this,  a  nonlinear  objective 

function  is  set  up  which  acts  as  a  rating  system  among  the 

feasible  candidates  and  selects  by  sons  criteria  a  "best" 

one.  The  criteria  will  be  to  find  a  top  of  spine  location  as 

close  as  possible  to  the  top  of  spine  position  at  the  beginning 

of  the  task  (denoted  as  (x  ,y  ,z  )). 

o  o  o 

(11)  Objective  Function  1:  (x  -  yQ)2  +  (y  -  y^)2  +  (z  -  zo)2 

This  method  has  the  advantages  of  ease  in  formulation,  and 
computational  rapidity  using  a  computer. 

In  general,  two  problems  are  solved,  each  corresponding  to  a 
protion  of  R3.  They  are: 

(12)  Problem  1:  Minimize  (objective  function  l) 

subject  to  Rl,  R2,  and  Pi,  i  =  1,2 

This  results  in  at  most  two  (true)  minimum  values  for  objective 
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function  l,  with  the  top  of  the  spine  located  in  at  most  two 
of  the  portions.  For  definiteness,  a  top  of  spine  location 
(x,y,z)  is  chosen  such  that  if  k^,  sire  the  minima  for 
problems  1  and  2,  let  k  =  MIN  (k^kg).  It  follows  that  (x,y,z), 
corresponding  to  this  minimum,  is  a  feasible  top  of  spine 
point. 


Corresponding  to  (12) 

(12a)  Minimize 

(x  -  xq)2  +  (y  -  yQ)2  +  (z  -  zq)2 


subject  to 

(la) 

r2  <  (x-x^2  +  (y-yx)2  +  (z-z^2  £  ^ 

(B2> 

(2a) 

r|<  (x-Xg)2  +  (y-y2)2  +  (z-z2)2  <  R^ 

(HI) 

(3a) 

r2  <  x^+y2  +  22<R2 

(4a) 

-R  sino;<  x  <  R  sin  a 

(PI) 

(5a) 

0<  y  <  R  sina 

(6a) 

R  cos  a  <  z  <  R 

(12b)  Minimize 


(x  -  xQ)2  +  (y  -  yQ)2  +  (z  -  zq)2 

subject  to 

(lb) 

ri  -  (x_x1)2  +  (y-yx)2  +  (z-^)2  < 

(»2) 

(2b) 

r2  <  (x-Xg)2  +  (y-y2)2  +  (z-z2)2  <  R§ 

($1) 

(3b) 

r2  <  x2  +  y2  +  z2 

(4b) 

(x-x^)2  +  (y-y^)2  +  (z-zg)2  <  B2 

(5b) 

-r  cos  <5  <  x  <  rcosfi 

(P3) 

(6b) 

0  <  y  <  rcosd 

(7b) 

r  sin  6  <  z  <  R  cos  a 

where  (x6,y6 

,Zg)  *  (A  sina  cosy  ,  A  sin  a;  sin  Y  ,  A 

0  <y<  n 

0  <oi<  n/2 
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Figure  12.  Region  of  Mutual  Intersection 
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3.2.2  Task  Modification 

If  a  task  is  infeasible  then  there  is  no  mutually  intersecting 
region  for  the  top  of  the  spine.  The  task  is  then  modified 
by  redefining  a  new  pair  of  palm  locations  which  can  be  reached 
with  the  link  configuration  under  consideration.  This  allows 
for  subsequent  evaluation  of  the  modified  task  which  "closely" 
resembles  the  original  one.  Ideally,  a  task  should  be  modified 
so  that  each  change  of  palm  position  could  be  an  indicator  of 
required  control  change.  This  would  provide  a  useful  design 
tool.  However,  relocation  of  a  control  must  be  considered  in 
relation  to  all  other  controls  which  may  be  operated  simul¬ 
taneously  with  it.  Thus,  it  is  beyond  the  scope  of  this  analysis 
to  provide  a  true  indicator  of  design  changes  on  the  basis  of  a 
given  task. 

The  modification  of  the  task  provides  that  a  top  of  spine 
location  be  calculated  so  that  the  resulting  body  position 
is  as  "close"  as  possible  to  the  original  task  defined  control 
locations.  With  some  alterations,  the  method  used  in  deter¬ 
mining  task  feasibility  can  be  used  here,  provided  that  the  left 
and  right  palm  locations  are  redefined  appropriately  so  that  a 
non-empty  mutual  intersection  region  exists.  It  is  important 
to  emphasize  that  the  two  palm  control  locations  of  a  task 
have  equal  importance  and  one  must  have  both  palms  positioned 
as  specified  in  order  to  consider  a  task  accomplished. 

The  cases  where  a  mutual  intersection  region  does  not  exist  are 
of  two  types:  (a)  R1  and  R2  have  a  mutual  intersection; (b)  R1 
and  R2  are  disjoint. 
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Case  I: 


INTERSECTION  OF  HAND  REGIONS 


This  case  allows  for 

RinR2PlR3  «0 
R1  n  R2  ^0 

It  is  necessary  to  find  the  top  of  spine  location  within  R3 
closest  to  the  intersection  of  R1  and  R2  and  once  chosen, 
designate  new  palm  locations  in  the  direction  of  the  original 
ones  and  at  a  distance  R^  from  the  top  of  spine  location. 

This  will  allow  BOEMAN-I  to  perform  the  task  in  the  «»■* 
manner  as  any  original  task. 

Suppose  (x^.,y^,Zj.)  represents  any  point  in  the  mutual  inter¬ 
section  of  Regions  1  and  2.  Let  a  spherical  surface  of  radius 
r  be  coincident  with  the  under-surface  of  Region  3  (see  Figure 
13).  There  exists  some  (x^,y^,z^)  whose  distance  to  the  sphere 
of  radius  r  is  a  minimum.  This  shortest  distance  is  embodied 
in  the  normal  from  this  point  to  the  sphere  at  (x,  y,  z) . 

It  is  determined  if  (x,  y,  z)  is  on  Region  3  or  below  it. 

If  (x,  y,  z)  is  on  Region  3,  then  let  (x,  y,  z)  *  (x,  y,  z) . 

If  (x,  z)  is  below  Region  3,  let  (x,  y,  z)  =(r  cos 6  cost], 
r  cos6  sinn  »  r  sin 6  )  using  the  rim  portion  of  Region  3  as 
the  circle  of  closest  location  to  (x^y^z^)  where  6=  tan 
and  is  the  angle  formed  by  (x,  y,  z),  (0,0,0)  and  (x,  y,  z). 

Thus,  mathematically,  let 

(13)  (x5t1)2  +  (y^)2  +  (z^)2  =  r2 
(x5t2)2  +  (y5t2)2  +  (z5t2)2  =  r2 

where 
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HK< 


Figure  13.  Case  I:  R1  O  R2  7^0 


D162-10128-1 

38 


Min 


and  represent  factors,  both  between  ±1,  which 

determine  two  points  on  the  sphere  forming  a  line  that 
extends  to  meet  (x,-,y,-,z^)  and  Is  normal  to  the  sphere. 


Choosing  the  t^  for  which 

^(x5-x5t1)2  +  (y^y^)2 +  (z5“z5ti)2  '  +  ^y5‘y5t2^2+^*5‘*5t2^2] 


holds,  let 

(x,  y,  t)  «  Min  (x5ti,  y5t±,  *5t±). 

If  (x,  y,  z  )  satisfies  the  constraints 

-r  cos6<  x<  r  cos6 
0  s.y  <  r  cos6 
r  sind^  z  < r 

then  (x,y, z)  «  (x,y,z)  is  the  chosen  top  of  spine  location, 
otherwise, 

(14)  let  (x,y,z)  =  (r  cos 6  cost]  ,  r  cos  5  sin  tj  ,  r  sln6), 

-  tan"1  (  £  ) 
x 

Thus  the  problem  is  to  find  an  (x^,y^,z^)  for  which 

(15)  (x^-x)2  +  (y^-y)2  +  (z^-z)2  is  minimized 
subject  to 

r2  ^(x5_x2^  +  (y5-y2)2  +  (z5_z2^2  -  ^  (Region  l) 

rl1^x5-xl^2  +  (y5"yi>2  +  ^25“zl^2^^  (Region  2) 
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MIN  | 

.  (x2"R2^(xl_Rl^ 

<  x5  <  MAX 

[  (X2+R2)>(xi+R1^] 

MIN 

[(y2_R2^^yi*Rl^] 

<y5<  MAX 

[  (y2+R2),(y1+R1)l 

MIN 

oT 

I 

rH 

N 

> 

1 

CVJ 

1 _ 1 

|  <  z5<  MAX 

[  (zg+RgMz^)] 

Using  (l4),  one  generates  a  feasible  top  of  spine  location  and 
redefines  (*2^2*  z2^*  (palm  control  positions)  to  be 


|*i  =  x  +  A^-x)*  Si 

y±  e  y  +  (y^y)* 

*i  =  2  +  (z-j-z)*  s1 


((xi-x)  +(yi-y)  +(z1-z2)‘1 ) 


1/2 


CASE  II:  DISJOINT  HAND  REGIONS 
Case  II  assumes: 


R1HR2  =  0 

It  is  necessary  to  find  a  top  of  spine  location  closest  to  each 
of  the  original  palm  locations  and  redefine  the  palms  as  in 
Case  I. 


Because  both  of  the  locations  hare  equal  weight,  the  midpoint 
between  them  is  chosen  and  the  top  of  spine  is  to  be  positioned 
within  region  R3  and  closest  to  this  mid- point  (See  Figure  13  ). 


If 

(17)  (x  ,  y  ,  z  )  = 
'  '  nr  m  m 


x,-x  y,-y? 

(*o  +  ),  (y0  + 


Z-| 

),  (z2  +  “§— ^) 


then  for  Regions  1,  2  and  portion  PI  or  P2  of  Region  3, 
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(x  ] >  Yy 


Figure  14.  Case  II:  R1  O  R2  -  0 


D162-10128-1 

41 


using  the  method  of  3.2.1, 

(18)  Objective  Function  2:  (x-x  )2  +  (y-y  )2  +  (z-z  )2 

m  m  n 

and 

(19)  Problem  i * :  Minimize  (Objective  Function  2) 
subject  to 

Pi' ,i'  =  1,2 

Choosing  the  smallest  of  the  two  resulting  minima,  the  top 

of  spine  location  with  respect  to  (x  ,y  ,z  )  is  found  and 

m  m  m 

Pi'  are  as  defined  in  (9),  (lO). 

The  task  is  redefined,  using  the  top  of  spine  location  (x,y, 
found  as  in  (l6). 
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THE  MOTION  MODEL 

This  section  mathematically  describes  the  model  for  the  simula¬ 
tion  of  human  motion  by  a  geometric  stick- figure.  This  method 
generalizes  easily  to  handle  rigid  skin  volumes.  The  problem 
is  formulated  in  a  general  way  in  the  Section  3. 3.1  where 
vector  functions  and  differentiation  formulas  relevant  to  the 
model  are  presented.  In  Section*  3.3.2, and  3.3.3,  the  mathe¬ 
matical  formulation  of  the  model  is  more  specific.  The  human 
motion  problem  is  presented  as  an  optimization  problem  with 
non-linear  constraints. 
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*-*> 


NOTATION 


B  =  a  vector  function  on  the  vector  z  of 

variable  Euler  angles,  defining  a  stick-nan 
body  configuration  in  the  vector  x,  of 
position  coordinates  of  link-connecting 
points  on  the  stick-man,  x  =  B(z). 

CI  =  a  Euclidean  coordinate  system  whose  center  is 

1  at  P_  ;  if  I  /  0,  the  z-axis  of  this  system 

1  i 

has  the  direction  of  PT  -  P_  ,  and  if 

-Li  1i-l 

1=0,  the  axes  of  C_  define  the  cockpit 
1  i 

reference  axes. 

F  =  an  objective  function  for  optimization. 

F(z)  =  f(z)  +  h(z)  where  f  and  h  are 
defined  below. 

=  the  composition  of  F  with  a  vector  function 
G,  which  is  defined  below,  to  get  an  objective 
function  in  terms  of  a  vector  y,  F(y)  «  F  (G(y)). 

=  an  objective  function  for  optimization  defined 
directly  in  terms  of  the  stick-man  body 
configuration  position  vector  x  =  B(z). 

=  the  composition  of  f  with  B  to  get  a  function 

directly  in  terms  of  Euler  angles,  f(z)  =  F(B(z)). 

G  =  the  vector  function  that  removes  angular 

constraints  from  the  stick-man  problem  by 
expressing  z  as  a  function  of  an  unconstrained 
vector  y,  z  =  G(y) . 

h  *  an  objective  function  for  optimization  defined 

directly  in  terms  of  the  Euler  angle  vector  z. 
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a  monotone  increasing  sequence  of  non- negative 
integers  defining  a  sequential  array  of  link¬ 
connecting  points  on  the  stick-man  (e.g.,  the 
Joint  locations  from  the  bottom  of  the  spiae 
to  the  top  of  the  spine,  then  on  out  to  the 
tip  of  the  right  arm). 


the  distance  of  I  P_  -  P 

Xi  xi-l 

length  of  link  Ii 


i 


i.e.  the 


a  3  dimensional  vector  giving  the  coordinates 
of  position  vector  Ij  on  the  stick-man 
(l  5  i  l  m) . 


a  column  vector  representing  the  point 

in  theC_  -  system. 
i 


a  non-negative  integer  such  that  is  an 

angle  for  rotation  T_ 

Xq(K) 


an  arbitrary  point  in  the  cockpit. 


a  column  vector  representing  S  in  the  C_ 

1i 


system. 


a  3-space  rotation  matrix  giving  the 

orientation  of  system  C_  relative  to  C_ 

ii  1i-l 


the  translation  vector  ,  expressed  in 

0 

0 


system  C_  , 

l 


used  to  translate  between  CT 

Xi-1 


and  C_  . 
xi 
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a  vector  such  that  z  *  G(y). 

a  vector  composed  of  all  the  variable  Euler 
angles  In  the  stictanan,  arranged  in  sow 
order . 

a  component  of  z. 

the  Euler  angles  used  to  calculate  T_  . 

1 
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3.3.1 


Coordinate  Systems 

Until  now  it  has  been  convenient  to  refer  to  the  locations 
P^(i  =  i,.  .  3*0  of  the  link  -  connectors  on  the  stick-man 

as  Euclidean  vectors  in  a  coordinate  system  based  at  Pq; 
however ,  a  means  of  expressing  the  orientation  of  each  body 
segment  (e.g.,  forearm,  hand,  etc.)  is  needed.  For  example, 
the  right  shoulder  (joint  l4)  not  only  haul  location  P^,  but 
also  an  orientation  that  expresses  the  amount  by  which  the  humeral 
link  -  Link  16  -  is  oriented  in  space  relative  to  the  right 
clavicular  link  -  Link  l4. 


To  express  the  orientations  of  Joints  and,  in  fact,  all  link  - 
connecting  points,  local  coordimtes  are  used  for  the  points  P 
on  the  stick-man.  Orientations  then  become  3- space  rotations, 
which  are  used  in  affine  transformations  to  change  the  system 
relative  to  which  any  point  in  the  cockpit  is  expressed  am  a 
Euclidean  vector. 


For  exactness,  let  ^  be  a  sequence  of  non-negative 

integers  defined  as  the  indices  of  the  points  P  in  a 

i 

sequentially  ordered  subsystem  of  the  stick-man,  e.g.,  the 
sequence  0,1,2,10,12,14,16,18,20,22  (where  m  =  10)  corresponds 
to  the  thorax-right  arm  system  (Figure  2).  Since  the  cockpit 
reference  system  is  based  at  PQ,  always  choose  1^  =  0. 

Define  coordinate  systems  C_  (i  =  1,...,  m)  centered  at 

■** 

P_  ,  each  with  z-coordinates  (z)  -  axis  parallel  to  link  I  . 

Ai  1 

Then  for  two  adjacent  points  I  .  and  I  ,  systems  C_  and 

1  1  ii-l 

CT  are  centered  at  a  distance  L  apart,  with  the 

i  Ai 

orientation  of  C_  relative  to  C_  expressed  by  the  matrix: 
i  1i-l 


T  =  T(9  ,  0  ,  *  ) 

li  lt  i  1i 


K  )). 

xi 
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©T  ,  0T  ,  and  ip  are  "quasi-Euler"  angles,  to  be  inter- 
Ai  1i  1i 

preted  as  indicated  in  Figure  1J>.  The  entries  for  the  natrix 

T(«,0,  -  (TJK(e,  0,  *))  are 


T 

11 

= 

cos 

6 

cos  0 

eosV 

+ 

sin  0 

siny 

T 

21 

- 

cos 

6 

sin  ^ 

cosy 

- 

cos  0 

sln,y 

T31 

33 

-sin 

e 

cosV 

T12 

« 

cos 

6 

cos  0 

siny 

- 

sin^ 

cosy 

T22 

- 

cos 

e 

sin  0 

siny 

+ 

COS  0 

cosy 

CM 

Eh° 

S 

-sin 

e 

siny 

T 

13 

- 

sin 

e 

cos  0 

T 

23 

- 

sin 

e 

sin  0 

T 

33 

a 

cos 

e 

t 

where  V  =  0  -  ip  . 
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X 


Figure  75.  Coordinate  System  Transformation 
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Thus  if  S  is  a  point  in  the  cockpit  expressed  in  system 
as  a  coluan-vector  S  ,  S_  is  found  By  Beans  of 

4  i-1 

the  affine  transformation  A  fire*  by 

*1 

ST  "  K  ST  =  T  (S-  +  t  )  vhere 

4-1  4  4  4  4  4 


0 


To  fora  S  ,  two  affine  transformations  are  applied  in 

4-2 

succession  to  get 

SI  m\  Ar  Si 

4-2  4-i  4  Ii 

"  Tt  (Tt  (s  +  t  )  +  t  ) 

4-i  4  4  4  4-i 


Ti  Ti  si  +Ti  Tr  4  +tt  4 

4-i  4  4  4-i  4  4  4-i  4-i 


Proceeding  in  this  way,  when  j  <  i,  then 

%  "THa  % 

T  s  +  / 

Ir  11 


K=J+1 


T  t 

4  4-q 


Note  that  if  S  is  P  ,  (the  coordinate  center  of  C  in 
4  4  4,  4 


C  -  coordinate^,  then 

4 


PI  ,1  =  M  4 
V4  k-j+1  4 


i4j+l)  i-q 

+  £  IT 

q=0  K=J+1 


4  4 

4  4-c 
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nay  be  expressed 


Letting 


PI  »  I 
i  3 


K 


n  V  PI  •  I 
:=i  TC  i  i 


i-Cj+D  <  . 

E  Vi  *1,  >  -  in  particular 

q=0 


PT  ■  PT  * 

Xi  Ii,° 


i-2 

z 

q=0 


Mi'q 


i-q 


The  expression  for  S_  in  terms  of  S_  now  becomes 

J  Ti 

si,  =  **3+1  st  +  V  ,1  5  and  tor  i  -  1*  I1  -  0) 

j  1  i  j  J 


S  = 
o 


3 

M 2  Sii  +  pj  >  an  egression  giving  S  in  Co  -  coordinates, 


In  pej-ticular, 

pi  *  ”s  h  +  PI 

i  1i  ii-l 


It  will  be  seen  in  later  sections  that  expressions  for  finding 
gradients  of  vector  functions  of  the  rotation  angles  0  , 

Ii 

^1^*  and  vj  »  discussed  before,  are  desired  for  the  present 

stick-man  model.  In  setting  forth  these  expressions,  only  the 
orientation  angles  that  may  vary  during  a  movement  of  the 
stick-man  are  treated  as  variables;  some  of  the  points  I±  have 
0  degrees  of  freedom  of  orientation  (e.g.,  for  I  =  3  or  5 
(Figure  2);  these  are  link-connectors  but  not  joints),  others 
2  or  3  degrees  of  freedom.  In  some  joints,  such  as  18  (see 
Figure  2),  the  angle  4>  will  remain  constant.  For  this  reason, 
a  parameter  vector  z  is  defined  whose  components  are  the 
orientation  angles  in  the  stick-man  that  may  vary  during  a 
given  motion  or  task.  If  only  the  right  arm  subsystem  is  to 
move  during  a  task,  the  components  of  z  are  chosen  as  the 
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variable  angles  in  the  right-arm  joints;  here  the  defining 
sequence  is  1^  =  0,  =  1,  1^  =  2,  1^  =  10,  I  =  12, 

I6  =  14,  I?  =  16,  I8  =  18,  I9  =  20,  and  I1Q  =  22. 


System  Cq  is  the  cockpit  reference  system  and  stays  fixed  by 
definition.  Referring  to  Figure  2,  fixed  orientations are  defined 
for  systems  C^,  Cg,  and  (centered  at  P.^  Pg,  and  P 

respectively),  and  for  (P^2  is  treated  as  a  link -connect or 
but  not  a  joint  for  the  present) .  For  the  body  segments 
centered  at  links  12,  lo,  18,  20  and  22  respectively,  2,  3» 

2,  2,  and  1  degrees  of  freedom  of  orientation  are  allowed. 

With  one  exception,  whenever  2  degrees  of  freedom  are  allowed, 
orientation  angles  8  and  <t>  are  used.  (The  z-axi3  of  system 

C_  is  rotated  counter-clockwise  an  angle  8  degrees  in  a 

i  i 

plane  degrees  (counter-clockwise)  from  the  x-axis  of 

l 

system  C_  .)  The  exception  is  joint  16,  where  0  « 
l-l 

and  are  used  (i.e.,  the  C^g  z-axis  is  rotated  8^ 

degrees  from  the  z-axis,  always  in  a  fixed  plane  in  Cl6, 

and  resulting  xy- system  is  twisted  counter-clockwise  by  an 
angle  ^g).  T^1US  the  parameter  vector  for  this  example  is 

z  =  (  ff12»  \2>  \6>  \6>  \6»  \q»  B2o>  02O» 


Given  a  component  variable  ^  of  t,  q(k)  is  selected  from  the 
domain  of  the  sequence  j 1^  }  such  that  z^  is  a  variable 

in  transformation  T_  (i.e.,  z,.  represents  one  of  the 

q(K)  ^ 

degrees  of  freedom  for  link  1^^) .  z^  apperas  as  a 

variable  only  in  T_  ,  thus  for  any  matrix  nroduct 

Xq(K) 


22) ' 
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the  partial  derivative  with  respect  to  z^.  is 


q(K)+l  if  i  <  q(K)  *  J 


K 


0  matrix  otherwise 

where  differentiation  is  done  componentwise  on  the  matrix 
T_  and  any  product  of  the  for*  is  interpreted  as  the 

\(K)  ^ 

multiplicative  identity  matrix. 


It  is  convenient  to  reformulate  (7)  for  recursive  differentia¬ 
tion  of  (i  =  2,..., a)  as 


i-1 


a z 


K 


■4’1 


&  K 
matrix 


if  q(K)  <  i 
if  q(K)  -  1 

if  q(K)  >  i  . 


Using  (8),  the  partial  derivatives  of  the 

point  locations  P_  as  fundi  one  of  the  parameters  z  can  now  be 

Ii  K 

expressed  according  to  the  recursive  formula 
A P  si  N  p 

_  I1  =  <***2  t_  +  *  *1-1  as  i  goes  from  2 

d\  1  d 


K 


to  m.  Use  has  bean  madeafthe  fact  that  t_  is  a  constant 

Ai 

vector  and 

A  p 

°  o  -  0  for  all  K. 


K 
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Gradient  formulas  for  vector  functions  of  the  parameter  vector 
z,  e.g.  the  gradient  of  the  point-vector  on  the  stick-man, 
will  now  be  exhibited. 

Suppose  that  z  has  n  component  variables  zR  (K  *  l,...,n). 

For  each  point  PT  in  a  system  defined  by  a  sequence 
Ai 

{ I  °=1,  formula  (3)  (or  (6))  defines  a  mapping  (vector 

function)  from  z  to  P_  .  This  mapping  is  differentiable, 

ii 

with  partial  derivatives  given  by  (9).  By  constructing  the 

,  where  (PT  ) .  is  the  j-th  component  of 
1i  J 


of  Pj  (j  =  1,2,3;  i  =  2,...,m),  it  is  possible  to  define  a 
1 

differentiable  mapping  B  :  z  — ►  x  =  B(z).  In  this  instance 
x  is  written  as 

rv 


(number  of  movable  points  in  the  system),  e.g.  in  the 
of  the  movable  right-arm  system,  the  vectors  Pq,  P^, 
P^q  would  not  be  used  in  constructing  x. 


P  =  3  x 

example 
P2,  and 


B^z) 


B^z) 


Bp(z) 


where 


vector  x  = 


(Vj 
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Corresponding  to  the  mapping  B,  for  each  value  z  for  which 

x  =  B  (z)  is  defined,  an  nxp  (gradient)  matrix  B  (i)  is  obtained 

z 

according  to 


)  -  ).  Each  entry 


(the  derivative  of 


has  the  form 


the  J-th  component  (j  -  1,2,  or  3)  of  P  )  for  some  i  and  J 

*i 

corresponding  to  r,  hence  may  be  found  using  (9)* 

The  preceding  analysis  will  be  uaedia  section  3.3.3,  where  a 
gradient-projection  optimization  model  for  the  stick-mam  is 
given,  in  the  following  way:  If  f :  x — *f(x)  is  a  differenti¬ 
able  function  on  the  x-space  defined  above  (e.g*  an  objective 
function),  then  for  each  x  »  B (i)  a  function  f  is  defined  cm 

z-spaoe  by  f(z)  =  f(B(z));  if  the  gradient  of  f  is  f  , 

yv  x 

the  gradient  at  f  at  z  «  z  is  given  by 


tz  (z)  0  B2(z)  fx  (x)  -  B2  (z)  fx  (B  (z))j 


A 

here  f  and  f  are,  as  usual,  column  vectors, 
z  x  9 
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3.3.2 


Constraints 


As  mentioned  before,  the  problem  of  simulating  human  motion 
is  that  of  finding  the  trajectories  of  all  body  segments 
when  a  task  is  defined  by  giving  the  trajectories  of  the  body 
segments  at  the  extremities.  Solutions  to  the  problem  are 
restricted  by  placing  limitations  on  the  stick-man's  motion 
due  to  human  body  constraints  (angular  constraints)  and 
environmental  constraints.  Thus  the  constraints  on  the  motion 
problem  are  of  two  kinds:  Constraints  that  serve  to  define  a 
task,  and  constraints  that  result  from  restriction  on  the 
freedom  of  movement  of  a  human  pilot.  In  terms  of  the 
mathematics  involved,  the  task-definition  constraints  are 
equality  constraints  and  the  limiting  constraints  are 
inequalities. 

TASK  DEFINITION 

The  present  formulation  of  the  task  definition  constraints 

includes  position,  orientation,  and  directional  constraints, 

all  of  which  have  been  used  in  the  existing  model.  Position 

constraints  require  that  certain  terminal  points  on  the 

body  reach  specified  points  in  the  crewstation  environment. 

This  is  the  most  important  type  of  constraint  in  defining  a 

task,  hence  it  is  used  in  all  cases.  To  defire  such  a 

constraint  explicitly,  suppose  is  a  terminal  point  on  the 

m 

stickman's  body.  It  is  required  that  P  match  the  point  C  in  the 

cockpit.  (C  may  be  either  the  control  m  point  to  be  reached 

by  Pj  at  the  end  of  the  task-motion,  or  an  intermediate  point 
m 

along  a  specified  path  from  the  initial  value  of  P^  to  the 

m 

control  point.)  P  =  P^  (z)  is  a  function  of  the  vector  z 

m  m 

of  Euler  angles  describing  relative  body  segment  orientations 
and  hence,  expression  (3)  in  Section  3.3.1  yields  the  constraint 
equation: 

(12)  PT  (z)  -  C  =  0  in  terms  of  z.  The  condition  (1.2)  is 
±  ra 
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satisfied  if  and  only  if  the  cockpit  reference  coordinates 

of  Pj  and  C  are  the  same, 
m 

The  orientation  and  directional  constraints  are  presently 
applied  only  at  the  end  of  a  task-motion)  to  more  completely 
define  the  terminal  position  of  the  task. 

A  more  precise  definition  of  terminal  position  in  terms  of 
orientation  seems  necessary  because  the  pilot's  body 
configuration  may  depend  on  not  only  the  location  of  the 
controls  he  is  to  reach,  but  also  the  orientations  of  his 
hands  when  grasping  the  controls.  To  formulate  this, 
consider  the  matrix  product 

■5  \  m- 

k 

This  is  the  product  of  all  rotations  T..  needed  to  find  the 

K 

coordinates  of  a  point  on  terminal  body  segment  I  in  terms 

m 

of  the  cockpit  reference  system  based  at  P  for  the  link- system 

defined  in  the  sequence  j  X.  (  M 

*  k*l 

MI  (z)  completely  defines  the  orientation  of  body  segment  I 

in  system  Co.  By  specifying  Euler  angles  0,  0,  and  4> 
corresponding  to  the  required  orientation  of  segment  I 

II 

this  orientation  can  be  defined  by  evaluating  the  rotation  matrix 

T(  0,  0,  if  )  =  f 

Hence,  it  is  required  that 

(13)  M  -  *  =  o. 

In  practice,  it  is  helpful  to  limit  the  number  of  entries  of 

(z)  and  T  actually  considered  as  constraints  (see  Section 
3.3*3) •  Presently,  only  the  first  column  of  the 
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matrix  equation  (13)  is  used,  giving  a  3-space  vector  equation 
which  expresses  the  condition  that  the  image  of  the  unit  vector 

under  the  rotation  (z)  matches  its  image  under  the 
rotation  T. 

Finally,  a  directional  constraint  is  applied  to  the  problem 
at  the  terminal  task-motion  position  by  the  requirement  that 
the  stick-man  be  able  to  look  at  a  specified  control  point.  It 
is  required  that  the  line-of-sight  match  a  vector  from  the  eye 
midpoint  of  the  stick-man  to  the  vision  control  point.  In 
the  present  model,  we  take  point  to  be  the  eye  midpoint 
and  consider  an  "eyeball"  located  at  as  an  added  link  of 
unit  length  based  at  P^.  This  link  has  two  degrees  of 
freedom,  corresponding  to  the  angle  ©  and  direction  0  of  bend¬ 
ing  of  this  link  away  from  the  z-axis  of  system  (aligned 
with  the  top  of  the  head).  In  the  present  model,  the  label¬ 
ling  system  has  been  modified  somewhat  for  ease  of  computing; 
the  line-of-sight  vector  is  P^  -  P^,  utd  it  is  required  that  this 
be  aligned  with  c-p6>  where  C  is  the  vision  control  point. 

Oae  way  of  expressing  this  is 

(110  (P7(z)  -  P,(z))  -  c  "  P6(z)  =  0;  for 

7  6  Kr-  P6(z)| 

purposes  of  experimentation  with  various  methods,  however, 
the  present  model  requires  the  inner  product  of  the  two 
unit  vectors  in  (lit )  to  be  unity,  i.e., 

(15)  (P7(z)  -  P,(z))  .  C  ~  P6(z)  -1=0. 

7  6  |C  -  P6(z)| 

This  completes  the  present  list  of  equality  constraints  as 
long  as  inequality  constraints  are  not  removed  by  redefining 
them  as  equivalent  equality  constraints  (see  Angular 
Constraints).  For  a  discussion  of  the  results  of  this 


1 

0 

.0 
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formulation  of  task  -  definition  constraints  to  date,  see 
Section  3.3*d* 


ABGRJLAR  CONSTRAINTS 

According  to  available  data,  the  constraints  on  the 
variable  orientation  angles  In  a  given  Joint  1  hare 
the  following  fora: 


That  Is,  the  lower  bounds  and  \  are  fixed,  but  the 

**1 

lower  bound  on  ©1  depends  on  0  ;  the  sane  relationship 
holds  for  upper  bounds. 


In  terns  of  the  paraaeter  vector  z  and  corresponding  stick** 
nan  sub-  system  described  by  the  sequence  j 1^  |  i*t 


Jq(K)  link  t0T  2k  i®  *  paraneter,  as  before. 

Thus  zR  Is  one  of 


0T  iB  a  variable  then  a,,,  *  ,  and  hence  the 

q(K)  rq(K) 


corresponding  inequality  from  (l6)  is 


(XT)  xK(Vl)  s  zK  s  mk(zk+1). 

If  z  is  either  0  or  ip  or  if  is  a  constant 

K  MD  V*)  Jq(K) 
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the  Inequality  takes  the  form 

(IS)  A,,  s  zK  i  uK  , 

irith  the  tipper  and  lower  bounds  fixed. 

As  a  means  of  unconstraining  the  problem,  suppose  a  differentiable 
transformation,  G,  is  constructed  that  takes  any  vector  y  with 
real  components  y^  (i=l,..,n)  to  a  vector  z  with  components 
z±  (i  -  1,  ...,  n)  in  such  a  way  that  the  z^  are  bounded  as 
above.  Thus  G  is  a  differentiable  mapping  from  real  n- space 
(y- space)  into  the  restricted  z-space  defined  by  inequalities 
of  the  forma  (17)  and  (l8).  If  it  also  is  required  that  G  is  a 
mapping  onto  restricted  z-space,  e.g.,  if  (17)  and  (18)  are 
taken  as  strict  inequalities,  the  component  mappings  G^ 
of  G  might  hare  the  form 

(17) '  zk  =  Gjy) 

'  \  (<Wy))  + 

=  (Gltl(y))  -  A  (G1+1(y))]  (i  +  (xAPWfy 

to  remove  the  constraint  of  equation  (17)  and 

(18) *  z±  =  Gi(y) 

■  K  +  (m  ~  \  )  (£  +  (l/ir)  Arctan  y  ) 
zi  zi  zi  1 

to  remove  the  constraint  of  equation  (18).  Equations  (17) •  and 

(l8) '  produce  an  unconstrained  problem  in  y- parameter  space 

(all  of  n- space ) ,  and  the  components  of  the  transformation  G 

3  !)•••)  n) • 
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Now  consider  the  function  f  defined  by  f(z)  =  f(B(z)), 
as  in  Section  3.3.1.  Defining  a  function  F  on  y-spece 
by  F(y)  =  f(G(y))  =  f(B(G(y)))  yields  the  gradient  formula 

F  (y)  =  Gy(y)  fz(z);  using  f^Cz)  as  given  by  (ll), 

(19)  Fy(y)  =  Gy(y)  Bz(z)  fx(x) 

=  Gy(y)  Bz(G(y))  fx(B(G(y))) 

for  y  in  y- space,  where  G  (y)  is  the  nxn  matrix  with  elements 

y 


If  the  mapping  G  has  the  form  given  by  (17) '  and  (18) ' , 

G  (y)  is  almost  a  diagonal  matrix. 

y 


An  alternative  means  of  removing  the  angular  constraints 
as  inequalities  is  to  transform  them  into  equivalent 
equality  constraints.  The  expression  for  an  overall  equality 
constraint  is 


(20)  g(z)  =  [H(Xi-zi)*(zi-Xi)2  +  HU^f^MZj-*^)2] 
i=l 

where  H  is  the  function  given  by 


H(x) 


j  0  if  x  %  0 
)  1  if  x  >  0; 


again,  if  some  of  the  constraint  boundaries  vary  as  in 
(17),  the  resulting  terms  in  the  summation  in  (20)  have 
the  form 


(21)  H(Xi(zi+1)-z.)*(zi-Xi(zi+1))2  + 

H(Vtti(*i+i)),(V,li(*i+i))2  • 


=  0 
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So  far,  both  formulations  (17)'  and  (20)  have  been  tried. 

The  added  complexity  inherent  in  using  (18) '  or  (21)  for 
variable  constraint  boundaries  has  eliminated  their  use 
in  the  Phase  I  motion-model.  Formulation  (17)'  has  been 
emphasized  because  it  automatically  eliminates  the  constraints 
from  consideration;  if  formulation  (20)  must  be  used, 
thar*  1*  om  more  equality  constraint  to  consider  in 
addition  to  the  task  definition  constraints. 

The  type  of  transformation  presented  in  (17)'  has  been 
found  to  work  only  in  a  one-arm  synthesis  model  used 
for  experimental  work  on  the  overall  motion-model. 

This  transformation  did  not  work  on  the  full,  upper 
torso  model,  a  deficiency  which  is  probably  a  result  of 
using  the  Arctan  transformation.  An  alternative  is  to 
use  sin  y^  in  (17)'  in  place  of  (h  +  (1/rr)  Arctan  y^); 
this  has  been  tried  on  the  one-arm  model,  with  performance 
considerably  improved  over  that  of  the  Arctan  version. 

Since  the  formulation  in  (20)  has  been  tried  with  no 

p 

success  on  the  upper  torso  model,  (17)'  with  the  sin' 
transformation  will  be  tried  on  the  full  man  in  Phase  II. 

ENVIRONMENTAL  CONSTRAINTS 

The  work  spare  of  the  pilot  will  generally  be  described 
as  an  enclosed  convex  region  bounded  by  planes  and/or 
spherical  segments  with  the  exclusion  of  certain 
protrusions  of  similar  geometric  description  (e.g., 
control  panels,  control  stick,  copilot,  etc©).  The 
environmental  constraint  implied  by  such  a  workspace 
definition  is  that  no  part  of  the  pilot's  body  may 
travel  outside  of  the  workspace.  These  constraints  are 
not  considered  in  Phase  I.  For  future  project  phases,  it  is 
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expected  tkat  eack  cob rax  cematraimt  be4y  caa  be  re- 
preaeated  by  mo  greater  tkaa  10  baale  eleaaata  (plaae 
aad  apkerlcal  aegaaata),  aad  tkat  tkere  id.ll  be  me 
■ore  tkaa  36  ceaatralat  bedlam. 
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3.3.3 


The  Objective  Function  Approach 


This  section  describes  the  application  of  optimization  theory 
to  the  simulated  man-motion  modeling  problem,  and  the  results 
obtained  so  far. 

Recall  the  parameter  vector  z  of  variable  Euler  angles  used  to 
orient  the  body  segments  of  a  computer  -  simulated  man.  A 
complete  description  of  the  entire  body  configuration  in  any 
given  position  is  obtained  once  the  corresponding  value  of  z 
is  kncsm.  In  the  optimization  approach,  this  is  obtained  by 
minimizing  F(z)  subject  to  the  constraints  of  section  3.3.2, 
where  F  is  a  real-valued  function  of  z  which,  when  minimized, 
expresses  a  set  of  hypotheses  about  human  motion.  (For  example, 

F(z)  might  be  the  total  energy  involved  in  maintaining  the 
body  configuration  corresponding  to  z.)  Using  the  penalty 
function  method  (see  reference  4)  to  remove  the  equality  constraints 
from  the  problem  (after  the  inequalities  have  been  removed  by 
either  of  the  two  methods  discussed  in  section  3.3.2),  the 
objective  function  F  is  modified  by  adding  to  it  a  penalty 
function.  The  Davidon  variable  metric  method  of  minimization 
(see  references  1  and  2)  has  been  used  with  some  modifications  to 
solve  the  stick  man  modeling  problem.  The  modifications:  a 
Fibonacci  search  (see  reference  3)  for  recovery  from  cubic 
interpolation  failure,  and  a  unimodality  test  to  ensure  success 
in  one  dimensional  optimization,  are  both  presented  in  reference  5. 

Suppose  the  objective  function  is  given  as  the  sum  of  two  terms, 

one  defined  directly  in  terms  of  position  coordinates  P  of 

Xi 

the  link-connectors,  the  totality  of  which  make  up  the  vector 
x  as  in  section  3*3 «1»  and  the  other  defined  in  terms  of  the 
vector  z  of  orientation  angles  for  the  stick-man  sub-systems 
under  consideration.  Let  the  two  terms  be,  respectively, 
the  functions  f  and  h.  Using  the  transformations  discussed 
in  section  3*3.1  an  expression  for  x  in  terras  of  z,  x  =  B(z) 
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is  obtained.  B  is  a  differentiable  napping  and  has  a  gradient 
matrix  given  by 

Bz(z)  =  where 


Xj  being  a  component  of  one  of  the  point-locations  PT  expressed 

in  C  -  coordinates.  ^ 

o 

The  function  f  :  x— f(x)  together  with  B  defines  a  function 
f  :  z — f (b(z) )  on  z.  Then  the  objective  function  evaluated 
at  a  point  z  in  parameter- space  has  the  form  P(z)  *  f(B(z))+h(z). 

Differentiating  as  in  formula  (ll) ,  the  gradient  of  F  at  z  is 
given  by 

(22)  F  (z)  =  B  (z)  f  (B(i))  +  h  (S) 

z  Z  X  z 

In  the  case  that  inequality  constraints  are  removed  by  use  of  the 
mapping  G  defined  in  section  3*3.2,  the  objective  function  F  in 
y-space  (defined  by  F(y)  =  F(G(y))  has  gradient 

(23)  Fy(y)  =  Gy(y)  F^z) 

The  method  of  differentiating  equality  constraints  (task  definition 

constraints)  for  inclusion  in  the  penalty  function  vill  be 

indicated  by  an  example.  Suppose  it  is  specified  that  a  terminal 

point  Pj  is  to  reach  control  point  C  (a  vector  in  C  -  coordinates), 
m  ® 
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As  in  section  3*3.2,  this  terminal  constraint  may  be  written 


(24)  Pj  (z)  -  C  «  0. 
m 

Define  the  equality- constraint  rector  to  be  g(z)  ■  (z)  -  C, 

m 

a  rector  with  three  components*  Using  the  transformation  0 
which  removes  inequality  constraints,  we  obtain  the  gradient 
matrix  for  the  vector  function  g  on  y-space  defined  by 

g(y)  ■  g  (G(y)): 


(25)  gy(y)  -  Gy(y)  gz  (G(y));  g2(z)  has  the  vector 


z=z 


as  its  K-th  row. 


To  sunmarize,  the  solution  to  the  stick-man  task-motion  problem, 
as  outlined  in  this  section,  is  to  optimize  on  the  equality- 
inequality  system: 


minimize 

F(z)  =  f(B(z))  +  h(z) 

I  subject  to 

g(z)  «  0 

and 

*zi  -  zi  -  ^zi 
or 

Xz±  ^zi+l^  ~  zi  -  ^  ^zi+l^* 

Using  the  constraint  -  removal  transformation  G,  we  may 
choose  to  solve  the  equivalent  problem 
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minimize 


F(y)  *  f(B(G(y)))  +  h(G(y)) 

(27) 

I  subject  to 

g(G(y))  =  0, 

where  G(y)  is  given  by  equations  (l?)’  and  (l8)'.  At  present, 
the  Davidon  variable  -  metric  method  with  penalty  function 
is  being  tested  as  a  gradient-projection  method  for  solving  the 
stick-man  problem  as  formulated  in  (27) . 


PHASE  I  OBJECTIVE  FUNCTION 

So  far,  the  objective  function  F  that  has  been  tried  with  this 
model  is  based  in  part  on  the  simplifying  notion  that  gravito  - 
inertial  effects  can  be  Ignored  in  formulating  human  motion 
criteria.  This  assumption  means  that  the  term  f(B(z))  is 
missing  from  F(z).  Thus,  F  required  no  evaluation  of  the  position 
coordinates  PT  (z)  that  make  up  the  vector  x  =  B(z).  These 

i 

coordinates  are  needed  only  for  the  constraint  vector  g(z). 

The  underlying  assumption  for  the  Phase  I  objective  function  F 
is  that  a  human  body  prefers  to  stay  as  close  as  possible  to  fetal 
(or  a  "Dead  Man  Floating")  position  at  all  times.  Thus,  if 
values  z^(i  ”  l,...,n)  are  selected  to  be  the  fetal-position 

Euler  angles,  the  squared  deviations  ( A  *^)  ■  (s^  -  z ”)  of 

the  angles  from  fetal  position  values  are  minimized.  Originally, 

the  form  of  F(z)  was 

n 

(28)  F(z)  =  ^  wi(z  -  z°)2, 

where  the  w  are  preassigned  constants  >  0  that  express,  for 
each  i,  the  relative  importance  of  keeping  ( A  z^)  as  small  as 
possible.  For  example,  the  weight  w^  for  a  deviation  Az^  in 
one  of  the  spine-system  (such  as  C^)  angles  would  be  considerably 
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greater  than  the  weight  w^  for  a  deviation  A z^  in  an  angle  that 
moves  the  radial  link  of  one  of  the  arms,  since  the  entire  upper 
torso  certainly  has  a  greater  tendency  to  minimize  muscular 
interference  than  a  radial  link. 

The  objective  function  in  (28)  was  used  with  success  on  the  one- 

arm  model  and  provided  a  basis  for  further  developmental  work 

on  the  fetal  position  idea.  Recently,  another  formulation  has 

been  tried  with  greater  success.  Let  A  be  the  set  of  Euler 

angles  used  in  making  up  z,  excluding  the  angles  0_  that  appear 

i 

as  variables  in  some  joints,  i.e.,  joints  where  the  direction 
of  bend  of  link  1^  is  variable. 

The  new  objective  function  is  as  follows:  If  z  corresponds  to 

a  "twist"  angle  </'  or  a  "bend"  angle  ©_  for  which  0_  is 
i  1  Ii 

fixed,  then  form  the  term 

(29)  Fr  =  If  Zy.  is  an  angle  for  which  0T 

varies,  form  the  expression  FR+1  for  zR+1 
as  in  (29) ,  and  for  z^  form  the  term 

(30>  fk  =  Vm  <**  - 

The  objective  function  is  then  formed  as  the  sum 

(31)  F(z)  =  Fi.  Here  an  attempt  has  been  made  to  express 
the  possible  dependence  of  F^  on  F^,+^  when  is  a  variable 

"bend-direction”  angle  0^  ,  instead  of  letting  Zj,+^  affect 
F(z)  independently. 
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It  my  be  that  the  angles  (ZL  should  be  removed  from  consideration 

altogether  in  formulating  an  objective  function  such  as  (28) 

or  (31);  there  my  be  no  preferred  direction-of-bend  0  for  a 

1i 


fetal  position.  Recent  results  obtained  by  setting  w  ■  0  for 
K  4  A  (z^  is  an  angle  )  indicate  that  this  is  so. 


Also,  a  better  way  of  expressing  the  preferred  orientations  of 
some  body  segments,  such  as  the  lumbar  region,  my  be  needed. 

For  example,  it  my  be  that  the  lumbar  link  has  more  than  one 
preferred  position  in  this  formulation;  e.g.,  if  it  is  necessary 
for  the  simulated  man  to  move  his  torso  forward,  a  re-defined 
preferred  Euler  angle  should  be  used  to  allow  him  to  settle 
more  readily  into  a  "hunched-over"  position.  Results  with  movements 
requiring  torso  excursion  strongly  indicate  this. 


In  conclusion,  further  work  on  the  objective  function  formulation 
should  include: 

(1)  continued  research  on  improvements  in  the  present  idea, 

(2)  research  on  the  validity  of  the  simplifying  assumptions 
mde  in  ignoring  the  effects  of  gravity  and  body-segment 
momentum  on  the  problem  system,  and 

(3)  reformulation  of  and  developmental  work  on  the  objective 
function  F  as  a  result  of  such  research. 

The  research  on  objective  function  formulation  outlined  above  is 
in  addition  to  that  needed  in  other  areas  of  analysis,  such 
as  improvements  to  the  optimization  method  being  used  and 
investigation  of  the  possibility  of  introducing  time- 
dependence  directly  into  the  problem  formulation  by  the 
use  of  additional  transformations. 
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3.4 


INTERFERENCE  ANALYSIS 


When  a  task  la  complete,  it  is  necessary  to  determine  whether 
BGEMAN-I  was  able  to  see  the  control  for  which  he  reached  arid 
whether  his  joint-link  system,  during  the  path  of  motion, 
encountered  the  work  station  geometry.  Interference  is  the 
general  term  describing  these  situations.  The  first  case  is 
concerned  with  visual  interference  due  to  the  intersection  of 
BCEMAN's  line  of  sight  with  one  of  the  cockpit  planes  and  the 
second  case  with  physical  interference  due  to  the  intersection 
of  at  least  one  link  with  these  planes  during  the  task. 

The  procedures  of  this  section  give  the  work  station  designer 
an  indication  of  visual  accessibility  of  controls  as  well  as 
the  dimensional  fitness  (to  a  given  pilot's  anthropometric 
characteristics)  of  the  overall  design  during  a  task  sequence. 

When  visual  interference  has  occurred,  it  is  desirable  to 
attempt  a  modification  of  BUEMAN's  point  of  view  to  regain 
visual  contact  with  the  control.  The  use  of  a  corrected 
line  of  sight  is  an  indicator  of  the  visual  complexity  of 
the  work  station.  Ideally,  in  any  work  station,  visual 
interference  due  to  the  environment  should  not  occur.  In 
practice,  environmental  interference  is  to  be  minimized. 

The  rest  of  this  section  is  concerned  with  detection  of 
visual  interference,  the  correction  for  this  interference 
and  the  detection  of  physical  interference  (restricted  to 
BOEMAN-I  and  the  seat  back  planes).  Correction  of  physical 
interference  is  not  attempted  during  Phase  I  of  the  project. 
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<«  M  4>  >  X  X  <X 


NOTATION 


Minimum  distance  from  intersection  point  to  a  plane  edge 

i,i=«l,2  =  Coefficients  of  Centroid  -  intersection  point  line  and 

plane  edge,  respectively 


i  =  1,2,3 

J 


(Vyj’V 

(x,y,z) 


line  and  plane  coefficients 
^th 

j  vertex  of  a  cockpit  plane 
a  point  on  the  line  of  sight 


a  point  on  redefined  line  of  sight 


A 

X, 


X 


M 


*MP 


intersection  of  Centroid  -  intersection  point  line  with  an 
edge  of  the  cockpit  plane 


(x5»y5,z5) 


(xc»yc»zc) 


(Vy4’V 

*<XM  *  XA) 


eye  aiming  point 

Centroid  of  a  plane  P 

Midpoint  between  the  eyes 

(x6,y6,z6)  =  Midpoint  of 
line  of  sight 


nearest  point  to  original  intersection  point  in  the  plane 
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3.4,1 


Visual  Interference 

For  the  baseline  model ,  the  Cockpit  Geometry  is  assumed  to  be 
composed  of  bounded  planes  which  are  polygons  with  four,  five 
or  six  sides.  These  planes  describe  the  panels  and  solid 
surfaces  of  the  work  station. 

A  line  of  sight  is  defined  with  the  end  points  given  by  the 
eye  midpoint  and  the  task  defined  control  point,  on  which  the 
eye  is  to  "focus". 

It  is  required  to  determine  if  any  of  the  bounded  planes  inter¬ 
sect  with  the  line  of  sight  and  if  so,  whether  it  is  possible 
to  define  a  new  line  of  sight,  by  moving  BOEMAN's  eye  mid¬ 
point,  for  which  no  intersection  with  the  offending  plane  will 
occur. 

These  problems  are  the  detection  of,  and  the  correction  for, 
visual  interference. 

DETECTION  OF  VISUAL  INTERFERENCE 

It  is  necessary  to  determine  if  the  line  of  sight  bounded  by 
the  eye  midpoint.,  calculated  by  the  motion  model  and  the  eye 
aiming  point  specified  by  the  task,  intersects  any  cockpit 
plane. 

Let  =  (x^,  y^,  z^ )  =  the  midpoint  between  the  eyes. 

=  (xq,  y^,  z^)  =  the  eye  aiming  point 
X  *  (x,  y,  z)  =  a  point  on  the  line  of  sight 

Then  the  equation  of  the  line  between  these  points  is  given 

by 


*  ■  (1  '  V  *  b 
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I  x  -  xh  +  SX5 

(l)  y  =  (l-t3)  y^  +  t3  y5 

(  x  *  (l-t3)  +  t3  z5 

where  0  *  t3  <  1 

Because  it  is  not  possible  to  easily  describe  a  bounded 
plane  with  four  or  more  vertices  in  analogous  form  as  (l), 
the  following  procedure  is  used. 

Let  Vj  be  the  vertices  of  any  plane  in  the  cockpit. 
The  are  ordered  such  that  is  adjacent  to  V^+1. 
Let  V1,  Vg,  V3  be  representative  of  the  plane  P  with 
VJ  =  zj)  for  J  2»  3.  (See  Figure  16) 


Then  the  equation  of  an  infinite  plane*  in  which  P  is  embedded, 
is  given  by 


x  '  Vi  *  ‘2t2  + 


or 


*  -  tjXj  +  tjXg  *  (l-tj-tj)  x3 

(2)  1  y  =  t1y1  +  t2y?  .  (l-tj-tj)  y^ 


I 


2  =  t121  +  t222  +  (X-tj-tj)  2j 


with  no  restriction  on  t,  and  t_ 

1  2 


The  line  of  sight  and  the  infinite  plane  are  tested  for  an 
intersection  point.  If  there  are  none,  the  line  is  parallel 
to  the  infinite  plane;  If  there  are  two  or  more,  the  line 
is  embedded  within  the  plane.  Otherwise  a  single  unique 
solution  exists  to  (l),  (2). 

♦Olmsted,  John,  M.  H. ,  Solid  Analytic  Geometry,  Appleton- 
Century  Crofts,  Chapter  2. 
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setting  (l)  and  (2)  equal  to  each  other  rearrang¬ 
ing  yield* 


(vi  -v3)  tx  *  (Va  -  v3)  t2  -  -  XA)  t3  -  Xn  -  v3 

/  (*1  -  x3)  t1  +  U2  -  x3)  tg  +  (x4  -  x6)  t3  -  x4  -  x3 


(3) 

(x3  -  yj)  tx  +  (y2  - 

y3)  *2  + 

(y4  -  y5)  t. 

(21  -  z3>  *1  +  <z2  - 

z3)  + 

(z4  -  z5)  t. 

Solving  by  determinants. 

let 

(xl  -  x3> 

(x2  -  x3) 

(x4  -  x5) 

(4) 

D  = 

(yl  -  y3> 

(y2  -  y3) 

(y4  -  y6) 

<Z1  ‘  z3> 

(z2  •  Z3} 

(z4  "  z5> 

K  -  x3) 

(x2  -  X3> 

(X4  ’  V 

(5) 

A  = 

(y^  -  y3) 

(y2  -  y3) 

(y4  ’  y5} 

K  ‘  z3> 

(z2  -  z3} 

(z4  '  z5^ 

(Xx  -  x3) 

(x4  -  x3) 

(x4  -  x6) 

(6) 

B  = 

(y2  -  y3) 

(y^  -  y3) 

(y4  -  y6) 

(zl  '  V 

(z4  "  z3} 

(z4  ‘  z5} 

<*1  ‘  x3) 

(x2  -  *3) 

(lU  -  X3) 

(T) 

c  = 

(y-L  -  y3) 

(y2  -  y3) 

(y4  -  y3) 

(zi  -  z3> 

(z2  "  z3} 

(zi»  -  z3> 
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(«' 


(provided  that  D  /  0) 


Tf  D  -  0,  then  either  no  solution  or  many  solutions  exist.  To 
distinguish  between  these  cases,  the  midpoint,  of  the  line  of 
s i eh t.  is  substituted  into  the  equation  of  a  plane  (in  determi¬ 
nant  form).  If  it  is  satisfied,  the  entire  line  mus*  be 
embedded .  Tf  it  is  not  satisfied,  the  line  is  parallel  to  the 
pi  ane . 


Tims 


XHP  "  +^XA  ’  XKF  *  W->’  y6’ 


is  t.he  midpoint,  of  the  line  of  sight. 


A  theorem  in  solid  analytic  geometry  states  :  "A  necessary 


and  sufficient,  condition  for  four  points  to  be  co-planar 
(x.v.Zj),  i  -  is  that  the  determinant 


=  0." 


x  y  z0  1 


?  ? 

xn  y,  zn  1 


1  J1  1 


*  Olmsted,  John  M.  N.,  Solid  Analytic  Oeometrv.  Apploton- 

r.ent.ury-Proft.s ,  p.  7)1. 
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Tf  (10)  is  not  satisfied,  the  line  is  parallel  to  the  plane, 
and  there  is  no  visual  interference.  It  is  expected  that  a 
line  embedded  in  a  plane  would  not  ordinarily  occur.  If  it 
did,  an  infinitesimal  movement  of  the  eye  midpoint  would 
eliminate  the  interference,  so  that  this  case  need  not  be 
considered  in  further  detail. 

Suppose  then,  that  D  /  0.  If  interference  is  to  occur,  it  is 
necessary  that  t^  in  (8)  satisfy 

0  <  t^  <  1  as  in  (1) . 

(t^  =  1  implies  that  the  intersection  of  the  line  of  sight 
and  the  cockpit  plane  occurs  at  the  eye  aiming  point  which  is  a 
normal  occurrence.)  If  not,  there  is  no  interference. 
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If  so,  one  oust  determine  If  the  intersection  has  occurred  vlthln 


the  bounded  plane  or  outside  It.  roar  this  the  centroid,  X  ,  af 

c 

the  plane  Is  calculated  using  the  vertices  Y  of  the  plane 

J 


(n) 


[e  ’  (V  ye'  ‘c)  aad 


J-l 


or 


D 

n 

n 

E  xy 

*e'£ 

%-E 

j-i 

J-l 

j-i 

D 

n 

n 

To  determine  If  the  Intersection  point  Is  vlthln  the  bounded 
plane,  a  line  Is  drawn  from  the  Intersection  point  X  to  the 
centroid  of  the  plane  X^,  known  to  be  Inside  the  bounded  plane. 
If  this  line  segment  crosses  suqr  one  of  the  edges  of  the 
plAns  (Vy  Vj+1)  then  X  Is  outside  the  bounded  plane. 


If  none  of  the  edges  are  crossed,  X  must  be  vlthln  the  bounded 
plane. 

Thus  ue  have 

(12)  Centroid  Intersection  Line 

2  -  +  x<j-(i-s1),  OS  sxs  1 

A 

where  X  Is  a  point  on  this  line 

(13)  Plane  Edge 


1  ■  VS2  *  Vl'(1-S2> 


whe  re  ^  i . 


0SS2S1 

J  “  * 

V  "V 
n+1  vi 


>,n; 


1  poin+  on  an  p-dre  of  a  plane . 
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If  no  common  solution  X  of  (12)  and  (13)  exists  for  any  j, 
then  X  is  within  the  bounded  plane  and  interference  exists. 
Otherwise,  there  is  no  interference.  Figures  16  and  17 
describe  each  situation.  Thus  in  all  cases,  one  may  detect 
the  occurrence  of  visual  interference. 

CORRECTION  FOR  VISUAL  INTERFERENCE  . 

It  is  necessary  to  distinguish  between  two  kinds  of  correc¬ 
tion  techniques:  The  general  and  the  specific.  The  general 
technique  provides  a  method  of  relocating  the  line  of  sight 
given  any  work  station  configuration  and  under  any  circum¬ 
stances.  A  specific  method,  described  below,  is  geared  towards 
a  simple  cockpit  geometry  in  which  the  occurrences  of  inter¬ 
ference  during  a  task  sequence  will  be  rare.  Under  this 
assumption,  the  situation  where  a  single  plane  causes  inter¬ 
ference  is  examined.  Visually  complex  cockpits,  where  two  or 
more  planes  simultaneously  cause  interference,  require  the 
more  general  method  for  solution  and  are  beyond  the  scope 
of  the  Baseline  Man-Model  for  Hiase  I. 

An  analysis  of  the  multimission  simulator  shows  that  the 
resulting  planes  are  connected  to  each  other  at  their 
boundaries  and  that  sources  of  interference  occur  only  for 
the  collection  of  planes  forming  the  "head-up  display"  and 
the  control  stick  platform  surfaces.  In  each  case,  inter¬ 
ference  would  occur  with  only  the  top  or  front  planes  of 
the  surface.  Thus,  it  is  sufficient  to  examine  a  corrective 
measure  for  interference  with  a  single  plane. 
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BOUNDED  PLANE  P 


Figure  16  Occurrence  of  Visual  Interference 
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Figure  17.  Visual  Interference-Free  Situation 


D162-10128-1 

80 


D 162-1 01 28-1 

81 


If  the  line  of  sight  Is  to  he  deflected,  it  is  desired 
that  it  he  moved  to  a  point  on  the  edge  of  the  plane 
vhose  distance  from  the  point  ef  intersection  is  a 
minimum.  This  situation  is  described  in  Figure  18. 

To  accomplish  this  it  is  necessary  to  find  a  line 
perpendicular  to  an  edge  of  the  plane  vlth  one  and 
point  at  the  intersection  point  X  such  that  its 
distance  D  is  the  minimum  of  all  possible  perpendicular 
lines  to  every  side  of  the  plane. 

Consider  a  plane  P  with  intersection  point  X  as  in 
figure  19, 

Let  P  have  vertices  Vj  (J  ■  1,2,..., 7)  vith  =  Y^. 

In  particular,  let 


(16) 


vhere  I  denotes  the  intersection  point  on  the  plane  and 
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PLANE  P 


Figure  19.  Determination  of  Minimum  Distance  to  Edge  of  Bounded  Plane 
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(18)  Dj  “  + 


Since  a  *4  not  available  directly,  aaetker 

* 

solid  analytic  gecnetry  is  invoked,  stating  If  a  Is 
tko  angle  between  two  directions  (\^,  ^  )  and 

(  X^,  ^2,  v 2)  then 


.  2 

sin  a  = 


*1  V1 

2 

'i  h 

2 

\  “i 

+ 

+ 

^2  v2\ 

v2  *2 

X2  »*2 

Let 


2j  •  <* 

M2.  -  (7 


J+l"^^4^ 

%  ■  <*lV/<2J 

The!  using  (17)  directly* 
2 


(7i-7J+i) 


(Tj-t  i)  (Tj"7.+1) 


J+l 

2\l/2 


then  cheesing  the  nlnlnun  D  -  Min  (D  ) 


The  location  of  the  nearest  point  T  on  the  closest 
edge  of  the  plane  say  be  described  a is  follows: 

A 

d„ 


(19)  let  u  -  ^2 

A 


*  Op.  Cit.  , 


P.  21 
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Then  Y  can  be  given  by 


(20)  Y  -  ulTj  *  (1  -  ux)  VJ+1 

Finally,  the  redefined  line  of  sight  Is  given  In  terns 
of  Y  and  the  eye  aiming  point, 

as 

(21)  X*  =  f  Y  +  (1-f)  XA;  f  >1 

where  x'  is  any  point  on  the  redefined  line  of  sight 
In  front  of  the  Interfering  plane  (towards  BCBtUM). 

Physical  Interference 

Since  BOBNAJf-I  mores  in  a  rigid  environment  It  Is 
desirable  to  have  an  Indication  of  when  link  inter¬ 
sections  occur  with  the  surrounding  geometry. 

For  Phase  I,  attention  Is  restricted  to  the  most  probable 
planar  encounter  -  the  seat  back.  In  general,  Wn* 
positions  are  given  and  presumed  consistent  with  the 
geometric  configuration  (no  Intersection).  However, 
Intersection  with  the  seat  back  Is  possible  for  all 
Interior  Joints  and  links. 

Future  phases  will  be  more  concerned  with  the  general 
detection  and  correction  of  physical  Interference. 

DETECTION  OF  PHYSICAL  INTERFERENCE 

The  problem  of  detection  is  handled  similarly  to  that 
in  Section  B.^-.l.  The  line  of  sight  corresponds  to 
a  link  (line  between  tvo  Joints)  and  the  seat  back 
plane  corresponds  to  the  cockpit  planes. 

Each  of  the  links  are  tested  in  turn  during  each  step 
in  the  task  and  all  seat  back  interference  is  determined. 
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Let  (X^,  Le  the  k**  link  between  points  X^,  Xfc+1 

let  Vj  denote  the  vertices  of  the  sent  back,  then 

(22)  X(k)  -  t^  +  (l-t3)Xk+1  Ogt3?l 

k  *  1,  «  .  «  ,  1 

is  the  equation  of  the  link  line,  where  m  is  the  total 
nonber  of  links  under  consideration  and 

(23)  X  «  ♦  t2V2  +  (l-t^-tgjv^  (where  t^  and  1 2  are 

suitable  constants)  is  the  equation  for  the  seat  back 
with  vertices  V1,  Vg,  V^. 

Setting  (22)  and  (23)  equal  yield* 

0*)  (v1-v3)t1  +  (v2-v3)t2  *  (Wi)t3  -  \+1-v3 

Using  determinants  to  salve  for  tj,  t2,  and  t^  similar 
formulas  as  in  (4)  -  (18)  are  generated.  Then,  if 

0<t3<  1  is  not  true,  physical  interference  has  not 
occurred,  (t^  =  0  or 1,  imply  only  that  a  Joint  is 
coincident  with  the  seat  back. 

Using  (ll)  to  calculate  the  seat  back  centroid  Xg,  the 
centroid  intersection  line  is  generated,  along  with 
the  seat  back  plane  edge: 


(25) 

x  =  xs1  xc  (i-s1) 

1;  k-1 

(26) 

*  ■  vjV  Vi(1-S2> 

0?S2^1  J-l 

y  =v 

n+1  1 
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where  a  Is  the  number  of  links  under  consideration  and 
X  represents  a  common  point  on  both  lines.  If  no 
oossnon  solution  X  exists  for  a mj  J*  then  X  Is  within 
the  seat  back  plane  and  physical  Interference 
occurred. 

Each  link  is  tested  in  this  fashion  and  the  resulting 
physical  interference  is  indicated. 

In  contrast  to  the  visual  interference  section,  no 
correction  for  physical  Interference  is  attempted  in 
I%ase  I  since  a  correction  procedure  would  actively 
involve  the  optimization  model.  An  overall  procedure 
for  determining  and  correcting  interference  during 
the  motion  is  scheduled,  for  Ffease  H. 
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